CUQM-133 



Comparison theorems for the Dirac equation with 
spin-symmetric and pseudo-spin-symmetric 
interactions 

Richard L. Hall and Ozlem Ye§ilta§t 

Department of Mathematics and Statistics, Concordia University, 

1455 de Maisonneuve Boulevard West, Montreal, Quebec, Canada H3G 1M8. 

Email: rhall@mathstat.concordia.ca and yesiltas@gazi.edu.tr 

^Permanent address: Department of Physics, Faculty of Arts and Sciences, Gazi 
University, 06500 Ankara, Turkey. 

Abstract. A single Dirac particle is bound in d dimensions by vector V{r) 
and scalar S{r) central potentials. The spin-symmetric S = V and pseudo- 
spin-symmetric S = —V cases are studied and it is shown that if two such 
potentials are ordered y(^) < V^-^\ then corresponding discrete eigenvalues 
are all similarly ordered E^J < E^^J . This comparison theorem allows us to 
use envelope theory to generate spectral approximations with the aid of known 
exact solutions, such as those for Coulombic, harmonic-oscillator, and Kratzer 
potentials. The example of the log potential V{r) = vln{r) is presented. Since 
V{r) is a convex transformation of the soluble Coulomb potential, this leads to 
a compact analytical formula for lower-bounds to the discrete spectrum. The 
resulting ground-state lower-bound curve Ei^(v) is compared with an accurate 
graph found by direct numerical integration. 



PACS numbers: PACS Nos.: 03.65.Pm, 03.65.Ge, 03.65.-w 

Keywords: Dirac equation, comparison theorems, spin-symmetric, pseudo-spin- 
symmetric, envelope method, log potential. 



R. L. Hall & 6. Ye§ilta§ 



2 



1. Introduction 

Recently there has been much interest in exact solutions of the Dirac or Klein-Gordon 
equations with scalar and vector potentials of equal magnitude [1-6]. By the term 
'vector potential' we mean the time component V{r) of the energy- momentum four- 
vector; the scalar potential S{r) is a term added to the mass. The spin and pseudo-spin 
symmetries in nuclear physics [7, 8] , which have been observed in the hadron, are used 
to explain aspects of deformed nuclei. Spin symmetry occurs in the spectrum of a 
meson with one heavy quark and anti-nucleon bound in a nucleus [9]. Pseudo-spin 
symmetry occurs in the spectrum of certain nuclei [10]. Ginocchio showed that the 
Dirac Hamiltonian with scalar and vector harmonic oscillator potentials admits both 
spin-symmetry and U (3) symmetry when S = V ,\i also admits pseudo-spin symmetry 
and pseudo-[/'(3) symmetry when S = — V in the three-dimensional {d = 3) case [11]. 
Exact solutions of Dirac equation with a Coulomb-like tensor potential is discussed 
in [12], under spin and pseudo-spin conditions. 

For non-relativistic problems where the Hamiltonian is bounded below and its 
discrete spectrum can be characterized variationally, the derivation of a comparison 
theorem of the form 

1/(1) < T/(2) ^ £;(!) < £;(2) 

is almost immediate. The situation is not so clear a priori for relativistic problems 
where the corresponding energy operator is not bounded below. In spite of this 
difficulty, it has been shown [13-17] that relativistic comparison theorems are indeed 
possible, at least with respect to vector potentials, that is to say, the time component V 
of a four-vector. To our knowledge, no such results have been obtained for comparisons 
involving a scalar potential S (a variable term added to the mass). In the present paper 
we derive comparison theorems for the Dirac equation in the case of spin-symmetric 
problems S = V and for pseudo-spin-symmetric problems S = —V. 

We consider a single particle that is bound by an attractive central vector and 
scalar potentials, respectively V and S, in d > 1 spatial dimensions and obeys the 
Dirac equation. For central potentials in d dimensions the Dirac equation can be 
written [18] in natural units 7i = c = 1 as 

ij-=H^, where H = 2_^asPs + {m + S)^ + V, (1) 

m is the mass of the particle, V and S are the spherically symmetric vector and 
scalar potentials, and {a^} and /3 are Dirac matrices, which satisfy anti-commutation 
relations; the identity matrix is implied after the vector potential V. For stationary 
states, algebraic calculations in a suitable basis lead to a pair of first-order linear 
differential equations in two radial functions {il)i{r),'ip2{r)}, where r = ||r||. For d> 1, 
these functions vanish at r = 0, and, for bound states, they may be normalized by the 
relation 

oo 

(Vi, Vi) + (V'2,V2) = J (V?(r) + V2'(r-))dr- = 1. (2) 



We use inner products without the radial measure factor A'^~^'> because the factor 

(d-l) 

r 2 is already built in to each radial function. Thus the radial functions vanish at 
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r = and satisfy the coupled equations 

E^^ = {V + m + S)i^i + {-d + kd/r)i}2 (3) 
Ei}^ = {d + fcd/r)Vi + {V-m- 5)^2, (4) 

where fci = 0, fc^ = T(j + ^^), d > 1, r = ±1, and the symbol d represents the 
operator d/dr. We note that the variable r is sometimes written u, as, for example 
in the book by Messiah [19], and the radial functions are often written ipi = G and 
•^2 = F, as in the book by Greincr [20]. We shall assume that the potentials V and 
S are such that there are some discrete eigenvalues Ek^v and that Eqs.(3, 4) are the 
eigenequations for the corresponding radial eigenstates. Here u is the number of nodes 
in the radial wave function for a given kd- In this paper wc shall present the problem 
explicitly for the cases d > 1. Similar arguments go through for the case d = 1: in 
this case k\ = 0, the states can be classified as even or odd, and the normalization (2) 
becomes instead (tl;f{x) + V'K'^)) dx = 1. 

The principal concern of the present paper is with two special cases, namely 
spin-symmetric problems, for which S = V, and pscudo-spin-symmetric problems, for 
which S — —V. In each class of problems, therefore, there is jiist one potential function, 
V: it is with respect to this potential that we shall derive comparison theorems. We 
shall first treat the two cases separately, and then, for the purpose of determining the 
discrete spectrum, put them together in a single formulation. 



1.1. Spin- symmetric problems S = V 
In this case (3) and (4) become 

+ — Vi = (m + i;)V2 (5) 
r 

^'^-^^2 = im+2V-E)^i. (6) 
r 

By differentiation and substitution we obtain the following Schrodinger-like equation 

-Vi' + ( M^^i±ll + 2{E + m)v) ^1 = -(m^ - E')^;,. (7) 



Although the radial function i/ji, which must be because of (2), is not separately 
normalized, the eigenequation (7) does determine the discrete eigenvalue E. This will 
be illustrated explicitly in the next section where we shall discuss some exact solutions. 
However, we note that (7) alone is not sufficient to derive the comparison theorem; 
for this it is still necessary to use the original Dirac equations. 

1.2. Pseudo-spin-symmetric problems S = — V 
In this case (3) and (4) become 

^[ + -^Vi = {m + E- 2V)i)2 (8) 

^2-— V'2 = (m-S)V^i. (9) 
r 

Again, by differentiation and substitution we obtain a Schrodingcr equation, namely 
^ (^M^zil + 2{E - m)V^ ^2 = -(m^ - i?^)^2. (10) 
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1.3. Combined eigenequation 

By using the following parametrization 

K = skd = ST + '^-^^ , IJ. = sm, s = ±1, (11) 

we may write the eigenequation for both cases as 

+ (^^^ + 2(^ + m)V^) = -(m' - E^)i;, (12) 

where tp is assumed to be square intcgrable on [0,00), and E = E^,^ is a discrete 
eigenvalue corresponding to a radial eigenfunction with u = 0,1,2, .. . nodes. 



1.4. Principal result 

The principal result of this paper is the following theorem: 
Theorem 

ya)<y(2)^i;a)<i;(2). (13) 

We shall prove this theorem in section 3. In section 2 below we first derive some 
exact spectral formulae for harmonic-oscillator, (;ouloinbic, Kratzcr potentials, and 
log potentials. In section 4 we outline the envelope method which provides energy 
bounds whenever a comparison theorem is available and the given potential can be 
written as a transformation V{r) = g{h{r)) of a soluble potential h{r), where the 
transformation function g{h) has definite convexity. In section 5 we consider the log 
potential V(r) = v In(; ) which is a convex transformation of the Coulomb potential 
h{r) = — 1/r. We show that envelope theory can use the exact Coulomb solution 
and the comparison theorem to generate an implicit analytical formula expressing the 
dependence of each eigenvalue on the coupling parameter v. The lower-bound curve 
El{v) is compared to an accurate graph obtained by scaling and direct numerical 
integration. 



2. Some exact eigenvalues 



We can obtain exact solutions by comparing the general Dirac eigenequation (12) 
for the classes of problems considered with the following generic radial Schrodinger 
equation 

'L{L+l) 



+ vf{r) ]i} = ei) = F^L{v)i}, 



(14) 



where u > is the coupling parameter for an attractive central potential with 

shape /(r), L is a generalized angular momentum quantum number that is not 
necessarily integral, and Fvl{v) describes how the eigenvalue corresponding to a 
radial eigenfunction with u nodes depends on the coupling. From the relation 
k{k + 1) = L{L + 1) we extract the following formula for L: 

1 

2' 



1 


1 




d-2\ 1 






ST H 


' 2 j + 2 


'^+2 


^ 2 ^ 





±1. 



(15) 



We now present six illustrations. 
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2.1. Harmonic oscillator: V{r) = vr"^ 
Prom Eq.(14) in this case we have 

£ = F,l{v) = Pv'^, P = (4J/ + 2L + 3). (16) 

If we apply this to the eigenequation (12) for the corresponding Dirac problem, we 
obtain 

E'^ - Ij? = {2v{n + E)Y P, ij. = sm = ±m. (17) 

Thus we must have v{E + fi) > 0, and we conclude \E\ > m. Numerical values are 
easily obtained as the sohitions to 

(E'^ ~m^)\E-ii\^AP^\v\. (18) 

satisfying \E\ > m. 

2.2. The linear potential: V{r) = vr 

For the linear potential we have by Eq.(14) and a scaling argument that 

S = F,L{v)=Pvi, (19) 

where the Schrodinger eigenvalues P = -F',yL(l) for unit coupUng v = I must be 
determined numerically; these are all positive. For example, if d = 3, r = s = 1, and 
j = |, then from Eq.(15), L = 1 : thus P = Foi(l) is the bottom of the spectrum of 
—9^ + + r, that is to say, P ss 3.3612545. If we apply this to the eigenequation 
(12) for the corresponding Dirac problem, wc obtain 

E"^ - fi"^ = P {2v{iJ. + E))i , fi = sm^±m. (20) 

Thus we must have v{E + jj,) > 0, and we conclude \E\ > to. Numerical values for 
Dirac energy E are then given by solutions to 

{E^ -m^){E- ij,f =4v'^P^. (21) 

satisfying \E\ > m. 

2.3. Coulomb potential V{r) = —vjr. 

We proceed in a similar way as for the oscillator problem. In the Coulomb case we 
have 

5 = F.i(z;) = -^, P = v+\^L. (22) 

If we apply this to the eigenequation (12) for the corresponding Dirac problem, we 
obtain 

(23) 

It is immediately clear that E"^ < = rn^, that is to say, the discrete eigenvalues 
satisfy —to < E < m. Moreover, we conclude that for spin-symmetric problems 
{jj. = to), we have v > but for pseudo-spin-symmetric problems (/U = —to), the 
coupling must be negative, v <Q. We note for future reference that in both Coulomb 
cases 

vji > 0. (24) 
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For the pure Coulomb problem we have the following explicit spectral formula 

E = ,x) ^. (25) 

(1 + ^) 

For the Coulomb potential we now briefly discuss the behaviour of and V'25 as 
solutions of (7) and (10), near the origin. The solutions of the combined equation (12) 
are given 



^(r) = cr^+ie-VA'^-^'^r^L+i (^^^2 _ , (26) 

where c is a normalization constant, >C°(x) are Laguerre polynomials, and the 
parameter L, which satisfies L{L + 1) = k[k + 1), is given explicitly by Eq. (15). 
The small-r asymptotic behaviour of the wave function for the corresponding 
semirelativistic problem is discussed in Ref. [21]. 

2.4- Shifted Coulomb potential: V{r) = —v/r + c 

We now add a constant term c and re-solve, obtaining, first from Eq. (12), 
where P = v + L + 1. Thus we have for the shifted Coulomb problem 



^ \ ' / 2c 2(u + c) 

"^(1 + ^) (1 + ^) ^ (1 + ^) 

We note that this result is consistent with Eq. (52) of Ref. [22]. From Eq. (12) in 
this case, it is clear that, for the existence of discrete eigenvalues, we must have 
v{E + ^) > 0. Prom this inequality and Eq. (28) we conclude 

v{c + /i) > 0. (29) 

Thus, unlike for the Dirac equation with a Coulombic vector potential and a constant 
scalar potential, where the Coulomb coupling must not be too large (the Z < 137 rule 
for atoms), here the magnitude |w| of the Coulomb coupling may be chosen as large as 
we please. Moreover, unlike for any Schrodinger problem, where an added potential 
constant term c may have any desired value, here c is restricted by Eq. (29). 

2.5. The Kratzer potential: V{r) = a/r^ — v/r + c 

The Kratzer potential [23-29] comprises a shifted Coulomb potential with an added 
centrifugal term. The corresponding eigenvalues are given implicitly by the same 
fomula (27) that we derived above, but now the parameter P depends on E. The new 
effective L parameter is determined by 

k{k + 1) + 2a{E + n) = L{L + 1). 

Since P = v + 1 + L, we can obtain (for a ^ 0) a second expression for E given by 



E=l- 
2a 



M- (30) 
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Thus, by equating the formulae (28) and (30) we obtain a quartic equation whose 
solution yields the value of P which then, in turn, determines E via (28) or (30). We 
can also use P from (28) which is 



P = vJ ^ + ^ (31) 
to obtain an eigenvalue equation which is equivalent to (30): 



v{ii + E) = ^u+^ + ^{K+^y+2a{n + E)j x 

{2c{n + E)+fi^-E^)K (32) 

We note that this result is consistent with Eqs. (40) and (52) of Ref. [22]. If c = 0, 
then the limit to Coulomb coupling v = yields from (28) E ~ ^ = sm = ±m, 
whatever value we choose for the centrifugal parameter a. We note that this conclusion 
clearly contradicts the spectral claims made in §(5.3.2) of Ref. [22]; we find no discrete 
spectrum if the Coulomb coupling is zero, that is to say when w = 0. 

2.6. The log potential: V{r) = v ln(r) 

We now let the potential shape be f(r) = ln(r) in Eq. (14). If we let the eigenvalue 
with coupling v = 1 and given v and L be 

e(l) = e,L(l) = F,L(l), 

then the general eigenvalue with coupling u > is given by 

F,L{v)^e{l)v-UMt,). (33) 

We see this by the following scaling argument. The Schrodinger Hamiltonian has the 
form H = — A + ?;ln(r). By scaling the radial variable r — >■ m, where cr > is fixed, we 
merely describe in a different way, and we obtain the spectrally equivalent operator 
H — >■ + vhi{ra). Thus the eigenequation becomes 



\ ] A + t;ln(m) 



il) = e{v)ip, (34) 



equivalently, 

[-A + £7^wln(r)] = a^{e{v) -vln{a))tp. (35) 

By choosing the scale so that cr^u = 1, the eigenvalue on the right-hand side of (35) 
must equal e(l). This establishes Eq.(33). If we now apply this formula to the Dirac 
combined eigenequation (12), we find that the Dirac energy E = E^l is given by the 
following implicit formula: 

E = n + v [2e(l) - ln(2) - ln(i;(/x + E))] . (36) 

For later reference we note that the bottom of the spectrum, for arbitrary v (of 
the appropriate sign) requires the single Schrodinger eigenvalue e(l) = ^oi(l) « 
1.6411353. From Eq. (12) we see that it is always necessary that u = + E) > Q. 
This parameter has the upper bound u < ui, where u\ corresponds to £■ = 0: if 
u> Ui, E becomes complex. We have from Eq. (12) with E = Q iov this case: 

-m^ = Ui{2e^L - ln(2)) - m \n{ui). (37) 
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Table 1. The four spectral regfons for the log potential V(r) = t;ln(r). 
When E = 0, M = v{E + fi) = v/j, = ui. 





II = 


HI 


// = — /;/ 


v>0 


-m < E < 


ui/v — m 


m < E < m + ui/v 


V <0 


—ui/\v\ — m 


< E < -m 


— ui/|d| + m < E < m 



For the ground state with p. = m = 1, d = 3, and w > 0, we have k = L = 1, 
eoi ~ 1.6411353, and Ui = u « 14.28389. More generaUy, depending on the signs of v 
and /U = ±m, we may identify four spectral regions: these are summarized in Table 1. 



3. Comparison theorem 

We establish the comparison theorem in two steps, the first establishes a differential 
result, as we did for pure vector potentials in Refs. [15,16], and the second extends this 
to a general comparison theorem, as in Ref. [17]. We therefore begin by considering 
spin-symmetric or pseudo-spin-symmetric problems in which V{r,a) depends on a 
parameter a and it is supposed that dV/da > 0. We shall prove below that this 
assumption impUes E'{a) > for each discrete eigenvalue. If we now define the one- 
parameter family of potentials by 

V{r,a) = V'^^\r)+a(v^^Hr)-V^^Hr)), ae[0,l], (38) 

then y(^)(r) > V^^\r) implies that V{r, a) is monotone increasing in the parameter a. 
Thus E'{a) > implies E^^) = E{1) > E{0) = E^^\ which result proves the theorem. 
It remains to prove the monotonicity of E{a). 

3.1. Proof that dV/da > ^ E'{a) > 0. 

If the normalization integral (2) is differentiated partially with respect to a, and we 
denote d'il)/da = tpa for each wave- function component, we obtain the orthogonality 
relation 

(V'la,V'l) + (V'2a,V'2) =0. (39) 

Now we differentiate (3) and (4) with respect to a to obtain 

E'{a)^i + E{a)^ia = (K + ^a)V'i + iV + m + S)^ia (40) 
+ {-d + kd/r)i;2a, 

and 

E'{a)ij2 + E{a)i,2a = (K - ^a)V'2 + {V-m- S)i)2a (41) 

+ {d + kd/r)i)ia- 

The linear combination (40)V'i + (41)V'2 of these two equations may be written 

£'(a)[(Vl, Vl) + (V'2, ^-2)] = (V^l, {Sa + Va^l) (42) 
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where 

W = (V2, [d + kdMiJia) - [ipi, {d - ka/r)i>2a) (43) 

-((m + E(a) + 5-y)V2a,^2) 

+ ((m-E(a) + 5 + F)Via,^i)- 

Wc have again used d to denote the differential operator d = d/dr. The boundary 
conditions imply the antisymmetric relation 

{^P,d(^) = -{di,,(^). (44) 
If (44) is used in (43), it becomes 

W = (Via, [(m - E{a) + 5 + y)Vi + {-d + kd/r)i^2]) (45) 
+ (V2a, [-(m + E{a) + S- F)V2 - {-d - fcd/r)Vi]). 

From equations (3) and (4) it is clear that W = Thus the expression for E'{a) 
becomes in the spin-symmetric case S = V 

E'{a)[{i,ui>i) + {i>2M = 2 (^Vi, I^V'i) (46) 
and in the pseudo-spin-symmetric case S = —V 

i;'(a)[(V'i,V'i) + (V2,V2)] =2(^V2,|^V2) . (47) 

Thus the theorem is established since, in either case, if dV/da has a definite sign, then 
E'{a) necessarily has the same sign. □ 
The Coulomb problem of §2.2 provides an illustration if we let a = w. We have for 
the shifted Coulomb potential V{r) = —v/r + c that dV/dv = — 1/r < 0. Meanwhile, 
we obtain from the corresponding spectral formula (27) 

_ 4(^ + c)v 

^^""^--p^ii + vyp^r ^^^^ 

Since, from Eq. (29), + c)v > for the shifted Coulomb problem, we conclude 
E'{v) < 0, as predicted by the comparison theorem. 



4. Envelope theory 

Envelope theory [30-33] is based on a simple geometrical idea. If the potential V{r) 
may be written as a smooth transformation V{r) = g{h{r)) of a soluble potential h{r), 
and the transformation function g has definite convexity, then the tangents to g{h), 
all of the form bh{r) + c, lie entirely above or below V{r). If the eigenvalue problems 
satisfy a comparison theorem, then these shifted /i-potentials provide upper or lower 
spectral bounds. The envelope method consists of finding the best energy bound from 
one of these families. Suppose, for example, that g{h) is a convex function {g"{h) > 0), 
then in this case we obtain a family of lower potentials given by 

V{r) > FW(r) = b{t)h{r) + c{t), (49) 

where 

b{t) = g'{h{t)), c{t) = g{h{t)) - h{t)g'{h{t)), (50) 
and r = tis the point of contact between the tangent bh{r) + c and the potential V{r). 
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5. The log potential 

We now consider the log potential V{r) = i'ln(r), where is a coupling parameter: v 
is positive in the spin-symmetric case and negative in the pseudo-spin-symmetric case. 
For a soluble envelope basis we choose the Coulomb potential h{r) = —1/r. Thus for 
the log potential shape wc have 

/(r) = ln(r) = g{h{r)) = - ln(-/i(r)). (51) 

It follows that g'{h) = —1/h > and g"{h) = > 0, that is to say, g is 

monotone increasing and convex. Equation (50) then provides us with the b{t) and 
c{t) coefficients: b{t) = t and c{t) = 1 + ln{t). The potential-shape inequality in this 
case therefore becomes 

ln(r) = f(r) > f^'\r) = + (1 + ln(t)). (52) 

r 

In view of the comparison theorem and the exact eigenvalue formula (28) for the shifted 
Coulomb potential, we are now able to construct spectral-bounds. Since the direction 
of the full potential inequality, inchiding the coupling, now depends on the sign of the 
coupling V parameter: if i; > 0, we obtain a lower bound; if t; < 0, envelope theory 
generates an upper bound. The bound, once established, may then be optimized over 
t. For our illustration we take v > for the spin-symmetric case /U = m, and we find 



El = max 
t 



(l- ^) 2v{l + ln(t)) 



p2 



(53) 



It is clear from this equation that as v ^ 0, El ^ /x, and as v ^ oo, El ^ —/-i. By 
working with a new optimization parameter q = (vt/P)'^, we can re-write Eq(53) in 
the form 

^ 1^\ ^ 2v{l + ln{P/v)) +vln{q) 



EL = ma^ + ^"^^ . (54) 

9 L Vi + 9/ 1 + 9 J 

After some algebra, we find that the critical value of q satisfies q — v/{p, + E), and, 
moreover, it is possible to obtain the following implicit analytical formula for all the 
eigenvalues: 

EL=m + v[l + 2 ln(P) - In {v{m + El))] , (55) 

where P — 1 + v + L and L is given by Eq. (15) . In Fig. 1 we exhibit the energy curve 
El(v) for the spin-symmetric groimd state (P = 2) for = m = 1, and in dimension 
d = 3, along with a corresponding accurate curve E{v) obtained from Eq. (36). 



6. Conclusion 

Comparison theorems immediately generate spectral approximations. Thus spectral 
data from one problem yields spectral estimates for another. Until recently, this 
kind of reasoning was not expected to be valid for relativistic systems since their 
discrete spectra are not easily characterized variationally. It was then discovered that 
comparison theorems could be established without recourse to variational arguments. 
This has led to some comparison theorems for relativistic systems with different 
vector potentials. In this paper we have shown that spin-symmetric and pseudo- 
spin-symmetric Dirac problems also admit comparison theorems. This in turn allows 
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5 10 15 

V 



Figure 1. Ground-state spin-symmetric energy E{v) for the potential V{r) = 
v\n(r), with m = 1, d = 3, and K = L = 1. The figure shows the envelope lower 
bound Ei^ provided by Eq. (55) and accurate numerical values E given by Eq (36) 
with e(l) = 1.6411353. 

US to use envelope theory to estimate the spectra of such systems. For the log potential 
wc arc thus able to generate a spectral formulae which provide energy bounds for each 
discrete eigenvalue. 
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